Surface Temperature Assimilation in Land Surface Models 
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Abstract - This paper examines the utilization of 
surface temperature as a variable to be assimilated 
in offline land surface hydrological models. Com- 
parisons between the model computed and satel- 
lite observed surface temperatures have been car- 
ried out. The assimilation of surface temperature 
is carried out twice a day (corresponding to the AM 
and PM overpass of the NOAAIO) over the Red- 
Arkansas basin in the Southwestern United States 
(31°50 'jV— 36°iV, 94 o 30W-104 o 30'kF) for a period 
of one year (August 1987 to July 1988). The effect 
of assimilation is to reduce the difference between 
the surface soil moisture computed for the precip- 
itation and/or shortwave radiation perturbed case 
and the unperturbed case compared to no assimi- 
lation. 

BACKGROUND 

Land surface modeling has faced limitations in the 
past due to the lack of spatially distributed data on 
land surface characteristics as well as variables in 
water and energy budgets, namely surface temper- 
ature and soil moisture. Soil moisture is a crucial 
component of both the water and energy budget. 
The absence of spatially distributed observations of 
soil moisture makes it very difficult for hydrological 
model validation. Comparison of model stream- 
flows at the catchment outlet with the observed 
streamflow does not ensure energy and water bud- 
get validation. Satellite observed surface temper- 
atures satisfies our requirements of being spatially 
distributed and having connections to both the wa- 
ter and the energy budgets. Surface temperature 


has seen observed by NO A A polar orbiting satel- 
lites using the AVHRR ( Advanced Very High Reso- 
lution Radiometer) and TOVS (Tiros Operational 
Vertical Sounder) since 1978. The subject of as- 
similation of soil moisture data or assimilation of 
meteorological data in order to estimate soil mois- 
ture more accurately is relatively a new area of 
study (McLaughlin, 1995). Recent advances in in- 
verse methods (Entekhabi et. al. 1994, Lakshmi et. 
al. 1997) have demonstrated the use of microwave 
satellite data in estimating soil moisture. The as- 
similation of soil moisture from low-level atmo- 
spheric variables using amesoscale model (Bouttier 
et. al. 1993a, b) have shown that the assimilated 
soil moisture estimates help in the initialization of 
atmospheric models. Another class of methods use 
satellite estimates of surface skin temperature to 
adjust for the soil moisture (McNider et. al. 1994, 
Ottle et. al., 1994) and estimate with greater ac- 
curacy the surface fluxes and surface temperature. 
Van de Hurk et. AL (1997) carry out assimilation 
by nudging the forecast model evaporation frac- 
tion using the satellite data and hydrological model 
computed evaporative fraction. 

THEORY 

The land surface hydrology can be represented by 
a tw r o-layer model as shown in Figure 1 (Mahrt and 
Pan, 1984; Lakshmi et. al., 1997). The water bal- 
ance for the model can be written as 
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where 6\ and 0 2 are the volumetric soil moistures 
of the top layer (with thickness zi) and the bot- 
tom layer (with thickness z 2 ), P is the precipita- 
tion, E is the bare soil evaporation, R is the surface 
runoff, T is the transpiration, q\ y2 is the moisture 
flow from layer 1 to layer 2 and q 2 ,wt is the moisture 
flow from layer 2 to the water table. In this model, 
the the transpiration is assumed to occur from the 
bottom layer only. The moisture flow from layer 1 
to layer 2 ( q \ j2 ) and the flow from layer 2 to the 
water table (q 2>w t) are modeled using the Philips 
equation accounting for the gravity advection and 
the moisture gradient. The difference between the 
model computed and the new evapotranspiration 
flux ET 1 (that satisfies the satellite observed sur- 
face temperature) is given by 


ET' - ET = SET = -4 «tT?6T s - H r ST s - G x ST s 

( 2 ) 

where ST s = T' s - T, , the difference between the as- 
similated surface temperature and the model com- 
puted surface temperature. The partition of this 
difference in evapotranspiration SET into the dif- 
ference for bare soil evaporation SE and the vege- 
tation transpiration ST is given by, 


SE = SET 
ST = SET 
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The soil moisture of layer 1 and layer 2 has to be 
modified by S6\ and S6 2 respectively so that this 
new bare soil evaporation and vegetation transpi- 
ration hold good. 



Figure 1: Representation of land surface hydrology 


This study was carried out over an areal extent 
of 4.75° in latitude and 10.5° in longitude in the 
southwestern plains of United States. The me- 
teorological data were obtained from the surface 
airways data (17 stations on a hourly basis) from 
Earthlnfo's NCDC (National Climate Data Cen- 
ter) data product. The TOVS data of surface air 
temperature and surface vapor pressure serve as 
input to the land surface model and the surface 
skin temperature is assimilated. The NOAA 10 
satellite with nominal overpass time at the equa- 
tor of 730am/pm is used. The assimilation is car- 
ried out for a period of 1 year (August 1987 - July 
1988). The land surface hydrological model is ba- 
sically run in two modes, viz without assimilation 
of satellite surface temperature data and with the 
assimilation of satellite surface temperature for the 
values of input precipitation biased 20% higher and 
lower. 
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where At is the time step in our land surface model. 


DATA AND METHODS 


RESULTS 

Figure 2 shows the comparisons between the volu- 
metric upper layer soil moisture averaged over the 
Red-Arkansas grid box between August 1, 1987 and 
July 31, 1988 for the control case with and without 
assimilation (top panel) and the difference between 
the control with assimilation and rainfall input de- 
creased and increased by 20% with and without the 
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surface temperature assimilation. 
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Figure 2: Surface air temperature measured in 

FIFE versus derived from NOAA 10 

It can be seen from Figure 2 that the assimilation 
of surface temperature from the NOAA 10 TOVS 
data leads to reduction in the volumteric soil mois- 
ture of the upper layer caused by incorrect rainfall 
input. In fact, the assimilation brings the soil mois- 
ture closer to the control case for both - reduction 
in rainfall and increase in rainfall. The ability for 
the surface temperature assimilation to serve as a 
reset is most important for global hydrological data 
assimilation. 
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